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GeV Emission from neutron-rich internal shocks of some 
long Gamma- ray Bursts 
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ABSTRACT 

In the neutron-rich internal shocks model for gamma-ray bursts (GRBs), the Lorentz 
factors (LFs) of ion shells are variable, and so are the LFs of accompanying neutron 
shells. For slow neutron shells with a typical LF of approximate tens, the typical /3- 
decay radius is ~ 10 14 — 10 15 cm. As GRBs last long enough [T 90 > 14(1 + z) s], one 
earlier but slower ejected neutron shell will be swept successively by later ejected ion 
shells in the range ~ 10 13 — 10 15 cm, where slow neutrons have decayed significantly. 
Part of the thermal energy released in the interaction will be given to the electrons. 
These accelerated electrons will be mainly cooled by the prompt soft 7— rays and give 
rise to GeV emission. This kind of GeV emission is particularly important for some 
very long GRBs and is detectable for the upcoming satellite Gamma-Ray Large Area 
Space Telescope (GLAST). 
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1 INTRODUCTION 

As realized firstly by Derishev, Kocharovsky & Kocharovsky 
(1999a), the fireball of Gamma-ray Bursts (GRBs) may con- 
tain a significant neutron component in several progenitor 
scenarios. For the neutron-rich fireball, the acceleration will 
be modified (Bachall & Meszaros 2000, hereafter BMOO) and 
there could be many novel observational signatures: (1) If 
the neutron abundance is comparable to that of proton, the 
inelastic collision between differentially streaming protons 
and neutrons in the fireball will provide us observable 5- 
10 GeV neutrinos/photons (Derishev et al. 1999a; BMOO; 
Meszaros & Rees 2000; Koers & Giannios 2007); (2) The 
very early afterglow emission of GRBs would be modified 
significantly because of the /3— decay (n — > p+e~ +u c ) of the 
"fast" neutrons carried in the outflow (Derishev et al. 1999b; 
Pruet & Dalai 2002; Beloborodov 2003b; Fan, Zhang & Wei 
2005b); (3) In the neutron-rich internal shock model, for 
GRBs with a duration > 20(1 + z) s, the (3— decay products 
of the earlier "slow" neutrons with a typical Lorentz factor 
(LF) ~ tens will interact with the ion shell ejected at later 
times and give rise to detectable UV/optical flashes (Fan & 
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Wei 2004; Fan, Zhang & Wei 2005a); (4) As shown in Vla- 
hakis et al. (2003), for a magnetized neutron-rich outflow, 
the neutrons can decouple from the protons at a LF V in the 
tens, while protons, which are collimated to a narrower angle 
by the electromagnetic force, continue to be accelerated to a 
value of r in the hundreds. As a result, a two-component jet 
might be formed (Vlahakis et al. 2003; Peng et al. 20050. 
Jin et al. (2007) claimed that the X-ray flat segment of short 
GRB 051221 was attributed to such a scenario. (5) Razzaque 
& Meszaros (2006) calculated the synchrotron radiation of 
the electrons resulting in the f3— decay and found unique 
temporal and spectral signature. 

In this work, we show that in the neutron-rich (un- 
magnetized) internal shock model, besides the possible 
UV/optical flashes predicted in Fan & Wei (2004), there 
could be strong GeV emission. This high energy radiation 
component is a result of the inverse Compton cooling of 
the electrons due to the spacetime overlapping between the 
prompt soft 7-ray photon flow and the interaction region 
between the proton shell and decay products of the neutron 



1 In this case, the protons and neutrons are separated in the nar- 
row and wide jet components, respectively. The decay products 
of the neutrons have no interaction with the protons collimated 
in the narrow core. If this is true, the discussion in this work is 
invalid. 
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shell. Such a cooling process, suppressing the UV/optical 
emission significantly, has been taken into account by Fan 
et al. (2005a). These authors, however, did not calculate the 
possible GeV emission. 



2 GEV EMISSION FROM NEUTRON-RICH 
INTERNAL SHOCKS 

2.1 Neutron-rich internal shocks 

In the standard baryonic fireball model, the prompt emission 
of long GRBs are likely to be powered by the violent collision 
of shells with different LFs (Paczyriski & Xu 1994; Rees & 
Meszaros 1994 ). Until now, the practical LF distribution 
of shells in neutron- free internal shocks model is unknown, 
let alone the case of neutron-rich. To reach a high efficiency 
of radiation, significant difference of LFs between two shells 
with approximately the same mass (M) are favored (Piran 
1999). 

The unmagnetized outflow is controlled by the dimen- 
sionless entropy r\ = L/(Mc 2 ) at ro, where L being the total 
luminosity of the ejecta and ro being the radius of the cen- 
tral engine. The decoupling of protons and neutrons in a 
fireball occurs in the coasting or accelerating regime, de- 
pending on whether r\ is below or above the critical value 
Tfr ~ 2.2 x lO 2 !^ 4 ^ 74 ^ + £)/2]" 1/4 , f being the ratio of 
number density of neutrons to protons (BM00). The conven- 
tion Q x — Q/10 x has been adopted in this Letter, in units 
of cgs. If a slow shell is concerned with rj 3 < ij^, neutrons 
are coupled with protons until the end of the acceleration 
episode. However, for a fast shell 0, rjf > rj^, neutrons de- 
couple from protons, keeping the velocity at that time, while 
protons can be accelerated to higher velocity. As a result, 
r„,s = Fp. a = r/s for slow shells and F n ,/ < V p j for fast 
shells, where the subscripts p, n represent the proton /neu- 
tron component; /, s represent the fast/slow shells respec- 
tively (BM00). So far we have four kinds of shells, i.e., slow 
neutron/ion shells and fast neutron/ion shells. 

Similar to the standard neutron-free fireball, the fast 
ion shell will catch up with the slower but earlier one and 
forms a new one moving with a LF T m ~ a few hundred at 
a radius i?i n t ~ 10 13 cm, suppose that the ejection time-lag 
between these two shells is about 0.01s (hereafter the formed 
"new" shell is named as the "i-shell"). The /3-decay radius 
of a neutron shell (hereafter, the "n-shell") reads as 



Re 



2.6 x 10 13 P„ cm 



For T n ~ tens — hundreds, Rp ^> i?i n t, so the /3— decay of 
these fast neutrons won't influence the physical process tak- 
ing place at i?i nt significantly. 

However, the much later ejected i-shell will catch up 
the earlier slow n-shell at a radius J? ca t w 2T^ l a cST/ (1 + z), 
where 8T being the ejection time-lag between the earlier 
slow n-shell and the later i-shell and z being the redshift of 
GRB. As long as R cat is comparable with R/3, s , i.e., 



ST^ 14(l + 2)(r„, s /30)" 



(2) 



2 In some afterglow modeling, initial Lorentz factors of the fire- 
balls ~ y/VfVs ~ 200 — 400 were derived (e.g., Zhang et al. 2006; 
Molinari et al. 2007). For rj a ~ tens, we have ~ 1000. So the 
assumption rj s < rj^ < rjf is reasonable. 



the decay- product of earlier ejected but slower neutron shells 
will be swept orderly by the later ejected ion shells in a range 
~ -Rcat — R/3,s and interesting UV/optical and GeV emission 
are produced. 

In this Letter, following Fan & Wei (2004), we assume: 
(1) The LF of shells has bimodal distribution, T c j = rjf or 
r e j = rj s with equal probability, which is favored by its rela- 
tively high efficiency of energy conversion and high peak en- 
ergy of the internal shock emission (Guetta, Spada & Wax- 
man 2001). (2) £ = 1. As a result, the ratios between the 
mass of the fast neutrons, slow neutrons and i-shells are 
~ 1 : 1 : 2. (3) Proton and neutron shells move with LF 
T pJ ~ 1000, r nJ ~ 200 for fast shells and r p , s = T n , s ~ 30 
for slow shells. After the merger of a pair of fast / slow ion 
shells, the resulting i-shell moves with (|Piran 1999T I 



MpjYpj + M p , a Tp, a 



200. 



(4) The ejecta expands into low density interstellar medium 
(ISM), as mostly found in afterglow modeling (Panaitescu 
& Kumar 2002). The ISM has been swept by the early and 
fast ions as well as the decayed products of fast neutrons, so 
the decay products of slow neutrons move freely. 

Notice that in front of the i-shells ejected at a time 
14(l+z)(r n , s /30) -1 sec, there are hundreds of decaying n- 
shells. With assumptions made before, the whole process of 
each i-shell interacting with these decaying n-shells is rather 
similar. For convenience, in our following treatment, the dis- 
crete interaction of each i-shell with these decaying n-shells 
has been simplified as the situation of an i-shell sweeping a 
moving proton trail (the LF of which is r n! „) with a number 
density (Fan & Wei 2004) 



2irR 2 m p Ri3 r 



(3) 



continually, where m p is the rest mass of protons and 
M n , 3 oc exp(— R/R@ >s ) is the mass of the decaying/slow neu- 
tron shell. 



2.2 dynamics 

The dynamical evolution of the i-shell is governed by the 
energy conservation of the system and can be estimated as 
(Fan & Wei 2004): 



(1) «7 



dm 



77rel 



Mi 



+ (l-e)U 



(4) 



where 7 represents the Lorentz factor of the decreasing i- 
shell and 7 ro i ~ (r njS /7 + 7/r njS )/2 indicates the LF of 
the i-shell relative to the n-shell. In the adiabatic situation, 
the radiation efficiency factor e is assumed as zero. U = 
(7roi — l)mc 2 is the thermal energy in the comoving frame. 
The swept mass satisfies dm = 4nnm p T n ^R 2 ([3^, — /3r„ >s )dR. 
f3~/ and /3r„ jS are the corresponding velocity of 7 and r n , s . 

In the numerical calculation, we take the following pa- 
rameters: the mass of i-shell and n-shell Mi = 2M n ~ 
5.6 x 10 27 g, the initial LF of i-shell T; w 200, the LF of n-shell 
F„, s =30, and the width of n-shell is taken to be ~ 3x 10 s cm. 
e e and es, the fraction of shock energy given to the electrons 
and magnetic field, are taken as 0.3 and 0.01, respectively. 
The minimum LF of the shocked electrons is estimated as 
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which gives (see also Fan & Piran 2006) 
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Figure 1. The thick solid line represents the evolution of the ex- 
ternal inverse Compton parameter (^eic)- Here, the luminosity 
of the prompt 7-ray emission is taken as 10 51 crg s _1 . The rapid 
increase of Y EIC is due to the exponential decline of the mass 
of the decay product of the neutron shell. Note that the upper 
thick line is the result assuming a series of parameters as in sec- 
tion 2.2, while the lower one is for the same parameters except 
£fl ~ £c = 0.3. The solid line indicates the dynamical evolution of 
the i-shell (7), while the dotted and dashed-dottcd line represent 
the minimum LF of accelerated electrons (7m) and the cooling 
Lorentz factor (7c), respectively. 
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2 m 



1 m 



^(7rcl-l) + l, 



(5) 



where m e being the rest mass of electrons and p being the 
power-law index of the distribution of accelerated electrons. 
Our numerical results are shown in Fig |2.2l 

The electrons lose their energy via synchrotron radi- 
ation and inverse Compton scattering on the soft 7— ray 



photons with a luminosity L 7 



10 51 erg s" 1 . 



The in- 



verse Compton is usually in the Thompson regime because 
7m-Ep/7 ~ E p < m e c 2 , where E p ~ 200 keV is the typical 
peak energy of GRBs (Preece et al. 2000). So the inverse 
Compton parameter can be estimated by 



where U n = L 7 /(47rJ? 2 7 2 c) and Ub 
l)nm p (? . The numerical estimate of Y E 



(6) 

rel 

has been pre- 
sented in Fig l2.2l One can see that the cooling of elec- 
trons is dominated by the inverse Compton process. Now 
the seed photons are from a region separated from the scat- 
tering electrons and along the direction in which the ejecta 
moves. Following Fan & Piran (2008) , we call this kind of in- 
verse Compton scattering as the External inverse Compton 
(EIC). The corresponding synchrotron self-Compton param- 
eter Y ssc can be estimated as follows. Following Sari & Esin 
(2001), in the case of fast cooling (see below for the reason), 
we have 



1 ^asr: ~t~ 



(7) 



Y 

1 ssc 



-(1 + Y EIC ) + V(l + y EIC ) 2 +4e c / 



0.2(1 + Y E 



for (1 + Y mc 
for (1 + Y E1C 

for (i + y EIC 



I 2 » 4e c /£B, 
l 2 «4e c / £B ,(8) 
I 2 < 4e G /e B . 



From Fig |2~2l we see that (1 + Y E1C ) 2 > 4e c /e B in most 
cases, so the EIC emission always dominates the SSC emis- 
sion. 

The cooling Lorentz factor y c of the electrons can be 
estimated as (Piran 1999) 



7c 



6nm e c(l + z) 
{l + Y mG )cr Tl BH> 



(9) 



where t is the observer's timescale, <jt is the Thomson cross 
section, and B is the comoving downstream magnetic field 
B fa v /327re e 7 re i(7 ro i - l)nm p c 2 . 

As shown in Fig |2,2l -y c <C j m . The electrons are in the 
fast cooling regime. 



2.3 GeV emission 

The electrons accelerated in the interaction between i-shell 
and the decay product of n-shells will subsequently Compton 
scatter the prompt soft 7— ray photons, and boost them to 
an energy 



h^„ 



7m Ep 



2 GeV(7 m /100) 2 (£ , p /200 keV). 



(10) 



The total energy of the GeV emission can be estimated 
as follows. Note that in our simplest model, the energy re- 
leased in an i-shell interacting with the j3— decay products 
of a series of neutron shells is equivalent to that of the in- 
teraction between an ion-shell and all the decay products 
of a neutron shell. So the GeV emission is largely the same 
as that resulting in the internal shocks at Rp, s powered by 
two ion-shells. Following Piran (1999), we consider a colli- 
sion between these two shells with mass Mi and M s which 
are moving with different LFs: r, 3> T n<s S> 1. Here, the 
same parameters are taken as in section 2.2. The resulting 
bulk LF r m in an elastic collision is 



ML, + M n r n , s 
Mi/r, + M n /r n 



108. 



(11) 



The total internal energy (the difference of the kinetic ener- 
gies before and after the collision) is 

E thc = r i A/ i c 2 +r„, 3 M n c 2 -r m (Af i +Af n )c 2 = 3.1xl0 50 erg.(12) 

Then the total energy given to electrons: 

E e = e e E the = 9.3 x 10 49 erg(^|). (13) 

Since the electrons are in the fast cooling regime, the 
GeV radiation efficiency can be estimated as (for Y EIC ^ 1) 



Y K 

1 EIC 
1 + ^ETn E to 



E, 



T.M.c 2 + F„, s M„c 2 



(14) 



With the typical parameters adopted in this work, we get 
?7 GoV ~ 8.5%. This is comparable to the typical GRB effi- 
ciency r) ~ 20% (Guetta, Spada & Waxman 2001). We thus 
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3 DISCUSSION AND CONCLUSION 



GLAST sensitivity 
time-integration of 30s 



10 2 10 3 10 4 

photon energy (MeV) 



Figure 2. The solid line represents the time-integrated GeV 
emission spectrum resulting from prompt soft 7-ray emission up- 
scattered by electrons in the neutron-rich internal shocks model. 
GLAST sensitivity with corresponding exposure time is also la- 
belled as the dashed line (Galli & Guetta 2008). 



draw the conclusion that the total energy of our GeV radia- 
tion component is comparable to the total energy of the soft 
7— ray photons in some very long GRBs. 

The EIC photons may be absorbed due to pair 
production by photons with energy above E a ~ 
2( 7 m e c 2 ) 2 /(h^Eic) ~ 1 MeV7 2 /(hz/ E ic/5GeV). And the 
corresponding optical depth can be estimated as (e.g., Svens- 
son 1987) 

11<7T-?V>E 



720tt_R 2 



0.2i?f 5 2 L 7 , 5 i(5T( 



En 



where N > E a 



200keV 

/3 7 L-.ST 



h^EIC - 

5GeV 



is the total number of pho- 
tons of the prompt emission satisfying hi/ > E a , ST ~ 
R l 3 iS /(2-y 2 c) is the timescale of the late prompt emission 
overlapping with EIC emission and its power-law index 
f3~f ~ 1.2 has been used to get the analytical coefficient. The 
ultimate optical depth is ~ 0.2 and the spectrum around 5 
GeV would not suffer a significant absorption. 

Now we estimate the detectability of this GeV com- 
ponent by the upcoming satellite GLAST. The Large Area 
Telescope on board covers the energy range from 20 MeV to 
300 GeV. The effective area around 1 GeV is S e n ~ 10 4 cm 2 . 
So the expected number of the GeV photons is 



Am 



47rD 2 hi/ EIC 



E~, 



10 52 erg v 1 GeV 



(15) 



Usually at least 5 high energy photos are needed to claim 
a detection (e.g., Zhang & Meszaros 2001), so we con- 
clude that this component is detectable for some bright long 
GRBs. Please see Fig |2.3l for a numerical example, which 
shows vF„ oc v 1 / 2 below the peak as a result of the scatter- 
ing electrons in the fast cooling phase. 



The physical composition of GRB outflows is not clear yet. 
The outflows, in principle, could be Poynting-flux domi- 
nated or neutron-rich. People found some evidences for ei- 
ther model. For example, the modeling of optical flashes in 
some GRBs suggested the magnetization of the reverse shock 
regions and favored the hypothesis that these GRB outflows 
were at least weakly magnetized, while the delayed/prompt 
optical flares detected in GRB 041219a and some other af- 
terglow modeling got some indication evidences for neutron- 
rich outflows (see Zhang 2007 for a review). The physical 
composition of GRB outflows may be not universal and 
more effects are needed to get some "unique" signatures of 
the Poynting-flux dominated or neutron-rich ejecta. In this 
work, we calculate the possible GeV emission signature of 
the unmagnetized neutron-rich internal shocks. This high 
energy component, only emerging in some very long GRBs, 
is accompanying some UV/optical flares (Fan & Wei 2004) 
and is expected to be correlated with the prompt 7-ray emis- 
sion but to have a ~ 10(1 + z) s lag. This is very different 
from Li & Waxman (2008), in which the UV/optical/GeV 
flares are powered by a series of weaker and weaker internal 
shocks at a distance ~ 10 lj cm and are almost simultane- 
ous with the prompt 7 — rays without an obvious lag. The 
time delay of the EIC emission, caused by the anisotropic 
emission of the scattered photons in the rest frame of the 
shocked medium (see Fan & Piran 2008 for a recent review) , 
~ (1 + z)i?/3 jS /(2 7 2 c), is also ignorable in the case we con- 
sider. 

In some very long bright GRBs, the delayed GeV emis- 
sion discussed in this work is energetic enough to be de- 
tectable for the upcoming GLAST satellite (see section 
2.3 for details). It, however, might be outshined by the 
synchrotron-self Compton (SSC) radiation of the regular in- 
ternal shocks at a radius J?i nt ~ 10 13 cm. So only for the 
regular internal shocks having e e ~ e_e, the delayed GeV 
emission may be the dominant signal. As shown in Fig |2.2l 
for e e ~ es, Y mc 2> 1 at t ~ i?/3 iS (l + ,z)/2 7 2 c > Is. So most 
of the energy generated at a radius ~ Rp, s is lost into the 
EIC GeV emission. For Gamma-ray Burst Monitor (GBM; 
< 20 MeV) onboard GLAST, the field of view is all sky not 
occulted by the earth and that of the Large-area Telescope 
(LAT; 30 MeV ~ 300 GeV) is ~ 2.5 sr. So ~ 1/5 GRBs will 
be detected by GBM and LAT simultaneously. With a good 
quality of keV— 10 GeV spectrum, we may be able to check 
our model. 
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